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ABSTRACT 



An extension of the transmission curves described in a previous paper is given. 
The effects of the earth's curvature and magnetic field are considered, as well as 
absorption or reflection from lower ionosphere layers. A chart is included for 
rapid calculation of transmission curves. 
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I. INTRODUCTION 

In a recent number of this Journal ! the author outlined a graphic 
method of obtaining, from vertical-incidence ionosphere measure- 
ments, the limiting frequencies and virtual heights of reflection of 
radio waves incident obliquely upon the ionosphere. In that paper 
was developed a type of " transmission curve" which, when super- 
imposed on a curve of frequency against virtual height, measured at 
vertical incidence, gave directly the virtual height of reflection for a 
given frequency at a given distance. A logarithmic sec <£ trans- 
mission curve was also described by which, for a given distance, one 
curve could be used to determine the heights, to a fair approxima- 
tion, at any frequency. 

It is the purpose of the present paper to evaluate in more detail 
the effect of the earth's curvature upon the ray theory of ionosphere 
transmission, to indicate the modifications consequently necessitated 
in the application of the transmission curves, and to consider briefly 
how the earth's field may affect limiting frequencies and virtual 
heights. 

II. EFFECT OF THE EARTH'S CURVATURE 

The influence of the earth's magnetic field will be neglected in this 
section. In the next section the effect of the magnetic field on radio 

1 Newbern Smith. Extension of normal -incidence ionosphere measurements to oblique-incidence radio 
transmission. J. Research NBS 19, 89 (1937) RP1013. 
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transmission will be considered for the case of the plane earth, and it 
will be assumed that the results hold approximately for the curved 
earth, also. 

The elementary theory of propagation of electromagnetic waves in 
an ionized medium whose surfaces of equal ionic density are planes 



"^ \tONO$PH£R£ 

EARTH 



Figure 1. — Variation, with height, of angle a straight line makes with the normal to 

the earth's surface. 

shows that, by SnelFs law, the waves will penetrate the medium until 
the refractive index \x r is reduced to the value sin <f> u where <f>\ is the 
angle of incidence of the waves upon the ionosphere, i. e., the angle 
the wave normal makes with the vertical on entering the ionosphere. 
To obtain the corresponding relation for a curved ionosphere it is 
necessary to consider the variation in the angle <j> which a straight line 
makes with the vertical, or normal to the earth's surface, at various 
altitudes (z). 

The geometry of figure 1 leads to the relation 

Hn*'=an*(l+££|^). (1) 

Assuming SnelFs law to be valid for a ray traversing an infinitely 
thin layer of the ionosphere, of thickness dz, we can use this relation 
to obtain the differential equation: 

d(n' sin 4i) dz 



\x' sin <f> R+h-\-z 



Integrating this from the lower boundary of the ionosphere, where 
/*'=!, <t>—<t>i, and 2=0 (see fig. 2) up to the altitude z, we get 

, . sin <fo 

M sm fr-— — (2) 

l ^R+h 
as the form of Snell's law appropriate to the curved earth. 
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This means that the wave will penetrate the curved ionosphere 
until the refractive index /x' is reduced to the value /V, given by 



Mo 



sin 0i 



14 



R+h 



(3) 



where <f>i is the angle of incidence of the wave on the lower boundary 
of the ionosphere and z is the maximum height of penetration above 
this lower boundary. 
From the geometry of figure 1, we have also 



sin </>i= sin 4>'( 1+^r^V 



(la) 



If we put <£'=</>o and z—z v , the values which these quantities have at 
the top of the equivalent triangular path in figure 2, this becomes 



sin <fo=sin ft/ 1+JJT^\ 



(4) 



TRANSMUTE 




IONOSPHERE 



Figure 2. — Transmission through curved ionosphere. 

z,= height of equivalent triangular path; *o=true height of reflection, at level where //=sln <j>i ( 1-^77 ) » 
R= radius of earth; and D= distance of transmission. 

and so the refractive index at the level of reflection must be 



/xo'= sin 0o- 



( i+ to) 



(5) 



instead of jUo /= =sin <f> 0) as is the case for the plane earth. 

If, as is usually the case, the penetration z^^R+h, we may write 
this condition 



(5a) 



686 Journal of Research of the National Bureau of Standards [Vol. go 

The effect of the earth's curvature is thus to cause the radio wave 
to be reflected from a region of smaller ionization density than would 
be necessary were the earth flat. 

In considering the problem of radio transmission over a distance, 
it will often be convenient to speak of an "equivalent vertical-inci- 
dence frequency" for the given wave frequency and transmission path. 
This is the frequency which is reflected at vertical incidence at the 
same height and undergoes approximately the same absorption as does 
the actual wave at oblique incidence. It will be denoted by /, and 
the actual wave frequency will be denoted by/'. 

Since, for a frequency/', 

where /o=-*/ — —> 23 it follows that reflection of this frequency will take 
place at a level where 



o' = A /l 



■f ' 



Z' 2 ' 



where /o' is the value of / at the level of reflection, corresponding to 
an ionic density of N . So 



/ 



lN e 2 



irm 



and/ is the frequency for which fi = at this same level of reflection, 
i. e., the equivalent vertical-incidence frequency. For a wave fre- 
quency/' the equivalent vertical-incidence frequency/, in the absence 
of a magnetic field, is therefore 



/=/Vi-mo' 2 , (6) 

where p 9 ' is the refractive index for the frequency /' at the level of 
reflection. 

In the case of the plane earth, we had 

/x ' = sin 0!=sin <£ , 
which gives 

J sec <£o 

the well-known secant law. In the case of the curved earth, however, 
we must use the values of /V given by eq 3 or eq 5, and 



V1-mo' 2 = L ■*** / 1 _ s j n .J_*±M ( 7 a) 



/l-Mo' 2 = 1 1 sin2 tl 




Under practical conditions z is almost always less than 400 km, 
and values of z 9 greater than 600 km contribute nothing to trans- 



i W. H. Eccles, Proc. Roy. Soc. (London) [A] 87, 79 (1912). 
1 J. Larmor, Jahrb. drahtlosen Telegraphie 25, 141 (1925). 
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mission over appreciable distances, so that we can always consider 
z and z v <^R-\-h. Using this approximation, eq 7a becomes 



Vl- Mo '2 =cos fcJi + ^- tan 2 c^cos fa^ l- 2( ^ +f f ttitffa(7b) 



Figure 3 gives values of 



vr= 



Mo 



-plotted against tan fa for various 
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Figure 3. — Variation of~7==== with z and tan d> . 
VI— Mo 2 

values of z , assuming h=100 km (at the bottom of the E layer). 
These curves may be used to determine the relation of / to / ' for 
various values of z and fa. By substituting fa for (j> x and z — z v 
for z , similar curves may also be plotted to determine the relation 
of j to/ 7 for various values of z , z V) and fa. 

If, further, z <^: (R+h) cot 2 fa or (z v —Zq) <C(J?+A) cot 2 <£ we 
may write eq 7b approximately, 

0x U=cos fa\ 1-^pHan 2 fa (7c) 



Vl-Mo' 



2 =cos fa\ 1 



Zo 



R+h 



tan' 



a form which it is convenient to use in some discussions. This ap- 
proximation leads to results good to 1 percent or better for jE-layer 
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transmission, where z Q is less than 50 km, and for single-reflection 
jP-layer transmission over distances less than 1,500 km, or mul tire- 
flection jF-layer transmission where each reflection covers less than 
1,500 km. For transmission over distances greater than 1,500 km 
for each reflection it is necessary to use the more exact expression. 

The correction term, or frequency A/ which it is necessary to sub- 
tract from/7 sec $o in order to obtain/, is the amount by which the 
equivalent vertical-incidence frequency for the curved ionosphere 
differs from that for the flat ionosphere. 

This is 

a/=/[cos 0o~ vr 11 ^] > 

which is approximately 

A/-/' cos *o[l-^/l-2 g=J tan 2 O ] 

and, if 

Zv— 2o<<(#+^) cot 2 O , 

A/ . .f,( Zy— Zo \ sin 2 <fr Q 
nj ~ J \R+h ) cos 0o 

This correction term is zero at vertical incidence and for z —2 v ', i. e., 
reflection from a sharp boundary at the height z Q . For z^z v ' it in- 
creases rapidly with distance. 

The vertical incidence frequency f derived above may be used to 
plot transmission curves of z v f against /, exactly as in the paper 
referred to in footnote 1. If the assumption be made that z v is the 
same as the 2/ measured at the frequency /, then these transmission 
curves could be superimposed on the curve of frequency against virtual 
height, measured at vertical incidence (hereafter to be called the 
(2/, j0 curve), to give directly the z v for the oblique-incidence case. 
The actual relation of z v to 2/ is discussed later. 

Actually, it is more convenient to use the^sec <frf transmission 
curves, wherein z v is plotted against f'/sec O , and to refer the correc- 
tion to the (z/,f) curve. On the assumption that z v =z v ', A/ can be 
calculated, for each point on the (2/, f) curve, for each sec O trans- 
mission curve. This transmission curve can then be superimposed 
on the (2/, f) curve displaced toward the higher frequencies by the 
amount A /, and the oblique-incidence values of z v read off. 

The calculation of A/ for each distance from 2/ and 2 , without any 
assumption as to 2„, and for a O calculated for an equivalent triangular 
path of height 2/, is discussed below. This value of A / can be used 
in the same manner as was described in the preceding paragraph, and 
the sec O curves applied, O being calculated for an equivalent triangu- 
lar path of height 2/ instead of z v . 

If the (2//) curve and the sec O curves are plotted logarithmically, 
as in the paper referred to, 4 we may calculate a factor 1+A///, as 
follows : 

1 i A /_ cos ^0 

* Extension of normal-incidence ionosphere measurements to oblique-incidence radio transmission, J. Re- 
search NBS 19, 89 (1937) RP1013. 
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The frequencies on the (z/,f) curve can then be easily multiplied by 
this factor by adding the factor logarithmically to the curve, and the 
log sec <£o curves can be applied to the resulting corrected («/,/) curve. 
It should be noted that in this expression the <£ refers only to an 
equivalent triangular path of height z/, and not z v . The z v enters 
only into the expression for /*</. 

A typical (z v ',f) curve corrected in this manner for a given distance 
(2,000 km) is shown in figure 4. The virtual heights are lower on 
the corrected curve and the curve extends out to frequencies higher 
than the critical frequency for the ordinary ray. The correction to 
the curve is different for different distances, increasing with the 
distance. This means that the virtual heights of reflection are lower, 
for greater distances, than if the curvature of the earth were not 
considered. The limiting frequency of transmission, or maximum 
usable frequency, over considerable distances is correspondingly in- 
creased, the difference becoming as much as 20 percent at the greater 
distances, depending of 
course on the form of the £ *°°r 
(z v ',f) curve. 

When the (z/, j) curve 
has been corrected in the 
above manner it may be 
plotted logarithmically, as 
described in the previous pa- 
per mentioned above, and the 
logarithmic sec <t> curves may 
be applied to give directly the 
maximum usable frequencies 
and virtual heights over trans- 
mission pathsof givenlengths. 

Since the first draft of this paper was written there has come to 
the author's attention an excellent unpublished paper "Skip-Distance 
Analysis," by T. L. Eckersley and G. Millington, in the form of a 
contribution to the November 1937 London meeting of the Special 
Radio Wave Propagation Committee held in preparation for the Cairo 
Radio Conference. In this they undertake an analysis of radio trans- 
mission over a curved earth and obtain curves for determining the 
retardation of sky wave over ground wave at a distance, in terms of 
vertical-incidence measurements. They limit the analysis, however, 
to the case where z is very small compared with h. The analysis can 
thus apply only to i?-layer transmission, since for i^-layer transmis- 
sion, z must be measured from the lower boundary of the ionosphere 
in order to include the effect of retardation in the E layer. 

Their work, however, suggests a method of approximate analysis 
for larger values of z and a means for easily estimating the error 
involved in the assumption made above that the z v for a given fre- 
quency /' was equal to the vertical-incidence z/ for the equivalent 
vertical-incidence frequency/. 

The development given by Eckersley and Millington begins with 
the following relation, for the ray path shown in figure 2 : 



.FREQUENCY -Mc/s 

Figure 4. — (z v ' f f) curve corrected for a given 
distance (2fiQ® km) for the effect of the earth's 
curvature. 
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where ds is an element of the ray path making an angle with the 
vertical. Since, from eq 2 

1 sin 0, 
sin 4>=— r 



M 1 , Z 



R+h 
this becomes 

ds 



,__sin_0i C B 



V^M + hy 

They also used the relation 



(8) 



C B ds = f- dz C» 

J A li' Jo fl' COS <f> J„ 



dz 




(9) 



R + h; 

For the case where there is no magnetic field ^^^jl—y^j where 
/«=.*/ — } and so 



•M /*' Jo 



/'<fa 



/ 4n_42_4t2 S131 ^1 M (\\ 



Putting into this expression the wave frequency f in terms of the 
equivalent vertical-incidence frequency given above in eq 6. 



C B ds = _ 1 f" /<fe = 

)a v.' vi^^^Jo rr~j7_ Zsin^_\ 



(id 



At this point Eckersley and Millington assumed that they were dealing 
with a thin layer, and that in consequence (1) 0' was only slightly less 

than -^zri ~~n from eq 8, (2) — was only slightly greater than 

1 f Z0 idz 
, I J > from eq 11. This simplification is justified for 

Vl~~Mo Jo V/ "~Jo 

J5-layer transmission, but is unfortunately not justified for i^ 2 -l a y er 
transmission, since the lower boundary of the ionosphere must be 
taken at the beginning of the 25-layer in each case. We shall not, 
therefore, confine our discussion to this limited case. It must also 
be noted that the simplification introduced by assuming 2 <C<C 
(R+h) cot 2 4>i also is not valid except at comparatively short distances, 
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where 0!=45° or less. For practical purposes, it is, however, justi- 
fiable to assume z <C<C B+A. Using this relation then, we obtain 

Eckersley and MiJlington combined the approximate forms of eq 8 
and 11 and introduced the value of V* - Mo' 2 =cos <f> u to obtain 

0'= o i i 1 I i% We shall, however, combine eq 11 and 12 to 

obtain the more exact relation 



where 



sin</>i f 20 fdz A _ 2z \ 

"(i?+A)Vr=^Jo v/ 2 d-^)-/o 2 V fi+A/' u ' a; 



i 
i-^i 



sim 0! 



fr+gjbJ i+ ^ tan2 ^ 



(12b) 



1-Mo' 2 ,. 3), o 



If we let 



#+A tan *> 



1+5 = 



V— 0-4?)' 

For values of 2 considerably less than (B+A) cot 2 <fo B is usually a 
small number, and for z— >2 (near the level of reflection) it is also 
usually small. It may attain, however, relatively large values. 

In general, the evaluation of B involves a knowledge of the dis- 
tribution of ionic density, and a precise evaluation of the quantity 
0' involves graphic integration for each case considered. We shall 
for the present consider two principal cases, z ->0 and z Q = z v ', the 
virtual height at vertical incidence for the frequency jf^/'Vl - Mo' 2 , 
and discuss some cases of linear distributions of ionic density, where 

2/Zo=/o7/ 2 . 

(a) 2 <c(i?+/0 cot 2 fa. This case applies to Z?-layer transmission 
at any distance or to J^-layer transmission at short distances (fa 
small). 

Combining eq 6 and 7 we get 

/=/' cos <fr[l +^i tan 2 fcl (14b) 

and from eq 12b 

1-^1+2(^5) taa**,. 
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From eq 12c 



1 



+ 5 = 1 +(ra)^J tan2 * 
f 



if we assume that (1 — z/z )/0- — fp 2 /f 2 ) is nowhere so large as to make B 
comparable with unity. This is equivalent to assuming that/ 2 // 2 
does not approach 1 much more rapidly than does z/z , an assumption 
that is valid except quite close to a critical frequency. Transmissions 
involving equivalent vertical-incidence frequencies close to the critical 
frequency are of interest, however, only when the distance of trans- 
mission is short, and in this case <£i is small, so that B is a small number, 
anyhow. 
Thus 

tan <t>T z , . "1 C zo dz 



s+iL x ~ra tan20i ]jo"/ 



Jl 



H 



R+h 



tan 2 <j) : 



'HO 



R+h 



Now J — / ./2 ==:gp/? ^e virtual height measured at vertical inci- 




dence for the equivalent vertical-incidence frequency /. Since the 
terms involving z /R+h are small, we may write 



g/ tan < fr ; 
R+h 



where 



t 1 -^ tan2 +*(i-Q-C'] ( 14c ^ 



7 _ 1 f» dz J 1 Z Q \ 

p ,__2__ f« zdz 

° «/(«+ A) Jo " 



For a linear distribution of ionic density 

-2z 0) and (7=1, and 



i- £ =i-i 2 >z. 



Zo 



/ 2 



,_tan_^i P» rfg T, _ 2g ~| 

"8+aJo 7jtTiL *+*-!' 

V Zo 
V^STa)} 



z/ tan <fr 



(15) 
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For z vanishingly small C'-*0 and eq 14 reduces simply to 

(b) z ->z/. This case applies to reflection from a fairly sharp 
boundary at the level z , the refractive index being nearly unity up to 
nearly this level. An example of this would be reflection from the 
sporadic E region. The angular distance of the part of the ray path 
where ju departs appreciably from unity is small and so the ionosphere 
can be considered as essentially flat. This is the condition assumed 
by the author in the paper already referred to (RP1013). Here 
z v —z/ and 

, . sin 6' , 1fi v 

tan O = — -, (16) 



B+h 
In this case 



■1 — cos 



Vl-/V 2 =cos o (l -%+fr tan2 *>) ( 16a ) 

For z/—z vanishingly small compared with (B+h) cot 2 O , this 
reduces to the case for the plane ionosphere, where 



VI — mo ,2 =cos 0o 
near distrib 
l—fo 2 /f=l — z/z and eq 12a becomes 



(c) Zq—^z/, i. e., a linear distribution of ionic density. Here 



PX 1 



2z_\ 
>'~ 2 Jo A 



e>== _ sin*, "V B+hJ 



-A-- 



Putting in the value of I— A and integrating 

/ sin^l+^tan^ ^ (17) 

Now since we may consider z^B+h we can write, as in eq 7b 

Vl-~^ /2 ==cos WHdti tan 2 fo 
and thus eq 17 becomes 

,__ z/ tan0 1 f 2 z/ 1 ( . 

6 —B+JTl^SB+hj (15) 

just as in the case where z <^(B+h) cot 2 1# We must, however, 
wiite for the equivalent vertical-incidence frequency in this case 

/«/ cos 0^1 +§^- h tan 2 X (17a) 
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instead of the value 

/=/' cos 0i[l +^ tan 2 X ], (14b) 

which we could use when z <^(R-\-h) cot 2 fa. 

The three examples just treated give an idea of what happens on 
transmission through the ionosphere, and of the angular distance 
0' the wave travels in the ionosphere as a function of angle of incidence 
fa (or vertex angle of equivalent triangular path O ), of true height of 
reflection z Q and of virtual height measured at vertical incidence z/\ 
We must now consider the part of the path from the earth to the 
lower boundary of the ionosphere. If we consider the incident ray 
to traverse an angular distance (B—0') in going from the earth to the 
lower boundary of the ionosphere, the geometry of figure 2 tells us 
that 

. J sin (0-0') , lft v 

tan fa= t s (lo) 

5+1— cos(0— tf) 

For a given fa, then, we may solve this equation for (0— 0'), and, 
by adding this angle to the angle 0' already computed, we obtain the 
entire angular distance traversed by the wave from the ground to 
the point of reflection in the ionosphere. For h—lQQ km, as we are 
assuming, it is sufficiently accurate to replace sin (0— 0') by (0— 0') 

and 1— cos (0— 0') by 77(0— 0') 2 in this equation. We can thus 



solve for (6— 0'), obtaining 



0-0' = cot 0i- J cot 2 fc~p (18a) 

We can now write the expression for the total distance of trans- 
mission D in terms of fa (or O ) ^o and 2/ : 



D=2R(cot fa- J cot 2 fc-ffV*'' < 19 ) 

where 6' is computed as above, from eq 12a. 

The time T required for the sky wave to travel from the transmitter 
to the receiver is, if c= velocity of the wave in vacuum, 

r V^ sin (0-0') C'ds\ 
c\ sin 0i ~ t J A fij' 

If we put in the value of J ds/fx obtained on the basis of the above 

analysis, and replace sin(0— 0') by its value in terms of b we may 
express Tin terms of the quantities X (or O ), z Q and 2/, as was done 
with D. 

The relation between the height z v of the equivalent triangular path 
and the virtual height 2/ measured at vertical incidence for the 
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equivalent vertical-incidence frequency /=/' JflTj^ can now De 
written. From the geometry of figure 2, 



tan 0o= tan (0i — #0 



sin 0' 



1— cos 6' 



so that 



R+h 
z,= (.fi+A)[sin 0' cot (^-flO-l + cos 6'] 



(20) 



The relation between z 9t z/, z , fa, O , $', D,f, and/ may be sum- 
marized for the cases discussed here. 



(a) Zo«^(R+h) cot 2 0! 



/=/' cos 0^1 + -^ tan 2 0,1 



r _ tan ^ f z f 1 p% 

6 -Ti+hl l ~~R+h^ *'JJ 



cfe 



7 2 



r=- 



fl(fl-g') , 1 p» 

sin 0! cos <t> : 



fo I 2z 



1J0 



V 



Jl 



20 



L-^ian^ l-_g 



{ 



1— ^ 

2() 



0— 0' = cot 0i — J cot 2 r 



2ft 



D=2R6, 

z v =(R+h)[sm 0' cot (fa — 6') — 1+cos 0'], 



(<0 ^o=0, 



"-« 



7« 



2 

I 

7 2 



/=/' COS 0!, 

/>/ _ g P / tan 0! 
e ~ R+h ' 



T _2V R(B-e') z/ 1 
cL sin 0! ~*~cos faj 



l-e'=o,ot fa- J cot 2 4>i-^r> 

D=2Rd, and 

z = (#+/*,) [sin 0' cot (0i — 0O -1+cos 0'], 
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(b) linear gradient of ionic density ( z =-^ z v ' j 

f=f> cos fo^i + ^L^ tan 2 X , 

,_ z v ' tan i r i _ 2 g/ "I 
*~ R+h L 3B+AJ 

cL sin </>! • R+hj 
0— 0'=cot 0i— -W cot 2 0i — o» 

Z>=2#0, 

2,= (B+A)[sin 0' cot (fa— 6')— 1+cos 0'], 

(C) 2b=2/i 

/=/' cos 0o, 

, , sin 0' sin 

tan 0o= — t 



-COS0' ^t^+l-cos0 



R+h ' A — i? 

r== 2 (g+ft) sin 

~ c sin 0o 
D=2R$, and 

Z^—- Z$ 

Case (a') corresponds to the special case treated by Eckersley and 
Millington, except that, for simplicity of analytical computation, they 
took R+h=R in the expression for 0', which introduces an error of 
about 0.016A percent (h expressed in km) — not a serious error for 
h=100 km, as we have assumed. This error becomes appreciable, 
however, if it is attempted to extend this treatment, as they have 
done, to a height of 400 km or so. 

Martyn's equivalence theorem, 5 developed for a plane ionosphere, 
tells us that 

Zv ==: 2v ) 

" %v * %% 

C COS 0o C COS 0i 

Z?'=2z/ tan O =22,' tan b and 
f=f cos O =/ / cos 0i, 
where 

0!= angle of incidence of waves upon the ionosphere, 

0o= half the vertex angle of the equivalent triangular path, 

z v =height of equivalent triangular path, 

T'=time the wave spends in the ionosphere, and 

D'— horizontal distance the wave travels in the ionosphere. 

We see that the relation z v —zj is valid in case (c) but not in any 

2 z ' 

other case. The relation T= '—r- is valid in cases (a f ) and (6), but 

c cos 0i v ' . 

is not valid in general. The relation D'—2z v ' tan 0i holds approxi- 
mately only in case (a'). And finally,/^' cos O in case (c), /=/' 

I D. F. Martyn, Proo. Phys. Soc. (London) 47, 332 (1935). 
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cos fa in case (a 7 ), and/ is a complicated function in every other case. 
It is therefore concluded that the equivalence theorem, in the form 
given, cannot be applied to the curved-earth problem. 

Referring again to the summaries of cases (a'), (6), and (c), we may 
for a given value of zj take a set of values of z , for each of which can 
be calculated the variation of D with fa, in each of the three cases. 
We may also, for these values of z ' and z , calculate the variation of 
/'// with fa. By eliminating fa graphically, we can determine the 
variation of/7/ with D for a given z v ' and z , or, what is the same thing, 
the variation of /'// with z v ' for a given D and z Q . We can plot a 
family of transmission curves, with/ as abscissas and z v ' as ordinates, 
for each value of D and several values of z , corresponding to each of 
these cases. 

Such curves may be superimposed on the (z, ', f) curve obtained at 
vertical incidence, and the (z , /) curve deduced therefrom, to give 
directly the wave frequency/' corresponding to reflection at a given 
level Zq, which is characterized by a given z/. The value of z appro- 
priate to reflection at the given height determines which curve of 
the family for a given D is to be used, and the point of reflection is 
determined as the intersection of this transmission curve with the 
(«/,/) curve. 

The value of z v depends only on D and fa } and so lines of equal z 9 
may be plotted on the curve sheet for each D, so that z v , as well as 
/', may be read off directly. 

Because of the approximations made in case (&) and the assumption 
of a linear variation of ionic density with height, this case is of only 
special significance. It will be assumed, until further investigation 
determines more precisely the variation of conditions with z 0j that the 
curves vary smoothly between those calculated for z — and those 
calculated for z Q =z v '. 

A family of curves for each distance is rather cumbersome for 
rapid use. It is, as was said above, more convenient to use the log 
sec fa transmission curves, and apply a correction to the (z v ',f) curve 
by multiplying each vertical-incidence frequency by the factor 
1+A///", where 

1 , A/ = cos fa 



7 V1-/V 2 

This factor is obtained, for a given D, z ' and z Q , by determining cor- 
responding values of D and V* — Mo' 2 for arbitrary values of fa or fa. It 
is unity for z v ' = z and is quite easily obtained for z —0. For inter- 
mediate values of z it will be assumed that the factor 1+A/// varies 
in a manner similar to that determined from the relation eq 7b with 
z 9 =z v , l. e., 

cos fa _ r 2(z v '-z ) . 2 " 

Vi-W 4 V s+r~ tan *° ' 

but drawn through the values for z =0 and z =z v / determined in 
the more precise analysis, rather than those indicated on the assump- 
tion that z v =z v '. Cos fa is here calculated for an equivalent triangu- 
lar path of height z/. 

57763—38 9 
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Figure 5 gives the approximate factors by which / must be multi- 
plied to give /' /sec 4> for values of 2/ from 200 to 500 km and for 
distances up to 4,000 km. 

III. EFFECT OF THE EARTH'S MAGNETIC FIELD 

The presence of the earth's magnetic field introduces some compli- 
cations in the use of these transmission curves. These complications 
are often of minor importance compared with some of the unknown 
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Figure 5. — Approximate factors by which the equivalent vertical-incidence fre- 
quency f must be multiplied to give the frequency f/sec <f>o which is to be used in 
conjunction with the log sec <f>o curves. 



factors (e. g., the geographic uniformity of the ionosphere over the 
transmission path), especially over long distances, but the effect of 
the earth's field must at times be taken into account. The anisot- 
ropy of the ionosphere due to this field causes the effect of the field 
on radio transmission to vary with the length, direction, and geo- 
graphic location of the transmission path. 

One effect of the field is to cause the received signal to be split in 
general into two main components, the one with the lower maximum 
usable frequency known as the "0" component and the other as 
the "x" component. The refractive index for a frequency/', in the 
presence of a magnetic field H, whose components along and trans- 
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verse to the direction of phase propagation are, respectively, i? L 
and H T is given 6 by 



n"=l 



ti 



f 



4'24 2 If //2/ 2 -]2 

/2_ J JT , / J JT I //2f 2 



where /n=-t/ > 

7rm 



W5 



^2^c Hl ' and 

2V= ionization density. 
Upper sign refers to o-component; lower 
to ^-component. 

The frequency of the wave whose x-component is returned from a 
given ionization density, at a given height, is different from the fre- 
quency of the wave whose o-component is returned from the same 
level. This frequency separation is in general a function both of fre- 
quency and distance and may of ten become negligible at great distances. 

For practical calculation, it may be assumed that only the field and 
direction of wave propagation in the region near the level of reflec- 
tion will appreciably affect the propagation of the wave. This assump- 
tion is probably better for the o- than for the x-component, and is, it 
must be emphasized, only a good approximation. 

With this limitation, therefore, the value of \x n may be calculated 
for a given transmission frequency and transmission path. If we put 
M /= =sin 0o and deduce z v from sin 4> as was done in the previous paper, 
we may plot transmission curves, of virtual height against equivalent 
vertical-incidence frequency, for the o- and x-components. When 
these curves are applied to the corrected (2/,/) curves, they may be 
expected to give reasonably good results. An example of this type 
of transmission curve is shown in figure 6. The one curve gives trans- 
mission conditions for the ^-component, and the other for the o-com- 
ponent. The frequency used is well over the gyrofrequency f H = 
eH/2wmc so that the ^-component is returned from a lower level than 
is the o-component and has a higher limiting frequency. 

It is not now justifiable to plot the transmission curves logarithmi- 
cally, since the form of the curves will vary with the transmission 
frequency /'. For practical purposes, however, a logarithmic curve 
may be used within a limited range of wave frequencies about the 
frequency j' for which the curve is plotted; a practical limit might be, 
say, within ±15 percent of this frequency. 

The logarithmic sec </> transmission curves may be used in estimat- 
ing the maximum usable frequency for each component over a given 
path by adding or subtracting the separation between the limiting 
frequencies for the two components, evaluated at that frequency and 
distance. This separation is in general a function only of the trans- 



• E. V. Appleton. J. Inst. Elec. Engrs. (London), 71, 642 (1932). 
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mission frequency and the quantity sec </> , and may be estimated within 
the limits of experimental error in most cases. 

Figure 7 gives the frequency to be added or subtracted from the 
maximum usable frequency given by the logarithmic sec O trans- 
mission curves for the cases H L =0, i? r = 0, and H T =H L . The H T =0 
case is uninteresting save for single-hop transmission over the mag- 
netic equator, close to the magnetic meridian. For transmission in the 
continental United States H L is much less than H T and, indeed, is 
negligible over east-west paths, so that such transmission we may con- 
sider as essentially transverse transmission. In this case, 



/ 2 =1- 



Jo 2 



-f'2 2 J JH 

J d /' 2 -/o 2 



where 5 = for the o-component and 1 for the z-component, and the 
logarithmic sec </> transmission curves give the correct maximum 





















S 


**c 










£ 














COM 


X 

PONF 




fKL 


O 


IFNT* 1 


/ 




X 
























/ 
/ 






o 

Id 
















f/i 


^-~ 




"nor* 


<AL-II 


ICIDE 


JCE 


X 
-J 

< 200 

1- 
k 
























CURV 


fc. 








S32523 


rf^ 






















> 
































6 7 

FREQUENCY 



MC/S 



Figure 6.- 



-Transmission curve, including effect of earth's magnetic field at level of 
reflection. 



Dotted line = secant-law curve (/'=/sec <£o) for the given distance D and wave frequency/'. Dashed line=» 
(*»',/) curve for the o-component, this being the component which, at vertical incidence is a measure of 
the ionization density N. The o-component of the frequency /' over the given distance is reflected at A, 
and the z-component at B. 

usable frequency for the o-component. For this reason the o-com- 
ponent lies along the sec 4> axis in the H L =0 case in figure 7. 

For a maximum usable f requency for the o-component below the 
gyrofrequency/ F = V/t 2 +/z. 2 the z-component always has a maximum 
usable frequency above j H > This is only important for i?-layer trans- 
mission and only in cases where the x-component is not too highly 
absorbed at frequencies near/#. This case must not be confused with 
the case of a transmission frequency/' less than f H , in which case the 
^-component is reflected from a level above the level where the o-com- 
ponent is reflected. 

Another effect of the anisotropy of the ionosphere due to the earth's 
field is to cause a difference in the directions of phase and energy pro- 
pagation in the medium. This results in the wave's being reflected, 
not at the level where the direction of phase propagation is horizontal, 
but where the direction of energy flow (group direction) is horizontal. 
This effect has not been considered in these curves, and is probably 
not important to the degree of accuracy to which these calculations 
are carried. 
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are independent of sec #o, save near sec 00= 1. 
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IV. BEHAVIOR OF THE WAVE BELOW THE POINT OF 

REFLECTION 

The development thus far has been given for any distribution of 
ionization in the ionosphere below the level of reflection. The effect 
on the reflection of the passage of the waves through the lower iono- 
sphere, and for that matter through the lower atmosphere itself, has 
been taken care of by using the vertical incidence virtual height in the 
calculation, since this effect enters into the measurement of this 
height. Nothing has been said thus far, however, about phenomena 
that happen below the level of reflection. 

An example of such a complication is given in figure 8. Here the 
transmission curve crosses both the i?and the F 2 virtual-height curves. 
Therefore, transmission may be expected by way of each layer. The 
angle of incidence is not the same for each layer, since the virtual 
heights of the layers are different. The wave reflected from the 
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Figure 8. — Transmission by way of both E and F 2 layers. 

Appoint of reflection on E layer; £=point of reflection on Fa layer; solid line = transmission curve; and 

dashed line= (z/,/) curve. 

higher layer has a smaller angle of incidence upon the lower layer than 
does the wave reflected from the lower layer. Under certain con- 
ditions, then, it can penetrate the lower layer, and the waves can reach 
the receiver by way of the higher as well as the lower layer. 

When the critical frequency of the lower layer is higher than the 
equivalent vertical-incidence frequency for the wave which would be 
returned from the upper layer, the wave will not penetrate the lower 
layer, but will reach the receiver by way of the lower layer only. 
Furthermore, if this equivalent vertical-incidence frequency is not 
below the critical frequency for the lower layer, but does lie in a region 
of appreciable absorption in this layer, the wave will be returned from 
the higher layer, but will be appreciably absorbed in so doing. Ex- 
perimental observations indicate the absorption of the actual wave to 
be roughly the same as the absorption of a wave of the equivalent 
vertical-incidence frequency measured at vertical incidence. 

It is possible to plot, on the sec <f> Q transmission curves, lines cor- 
responding to the values of sec <f>i (see fig. 2) for different heights. 
Such a set of "absorption lines" is shown in figure 9. When the 
transmission curve is superimposed on the (2/, f) curve the behavior 
of a wave below the reflection level may be estimated by the region 
of the (2/ , f) curve sheet through which the absorption line through 
the reflection point passes. If this line passes through a region of 
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absorption or cuts a lower layer the wave will be absorbed or will not 
penetrate through to the higher layer. 

These absorption lines are the lines sec 0!=const. for a plane earth 
and for short distances on a curved earth. They curve toward larger 
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FREQUENCY Mc/S 

-Absorption lines on logarithmic transmission curve. 



Dashed line=(z/, /) curve; dotted lines = absorption lines. E reflection takes place at point A and Ft 
reflection at point B. Fi reflection at Cis shielded by the E layer at D. 

values of sec $j for lower heights in the case of greater distances over a 
curved earth, and approach the transmission curve itself for great 
distances. An example of the use of the transmission curves and 
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Figure 10. — Virtual heights at vertical and oblique incidence. 

(a) (z/,/) curve showing E, Fi, and Fi layers. (6) oblique-incidence heights corresponding to (a). Dotted 
lines show parts of curves shielded by E layer. Note transmission of some frequencies by several layers. 
This shows how the i*Vlayer reflections are relatively unimportant for transmission, (c) (?.',/) curve 
showing absorption (X) above critical frequency for the E layer, (d) oblique-incidence heights corre- 
sponding to (c). Absorption here (X) completely blocks out some frequencies. All curves are for the 
o-component. 

absorption lines is shown in figure 10, which gives vertical-incidence 
curves and the corresponding frequency-height curves derived there- 
from for the transmitted wave over a distance. From these latter 
curves can be deduced, of course, the equivalent paths or group re- 
tardation of the waves over the transmission path. 
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V. ANGLE OF DEPARTURE AND ARRIVAL OF THE 

WAVES 

Disregarding any possible asymmetry of the wave trajectory due to 
the earth's magnetic field, the angles of departure of the waves from 
the transmitter and of arrival of the waves at the receiver are equal. 
The amount of energy which is radiated from the transmitter at a 
given angle from the horizontal depends, especially for low angles, 
on the design of the transmitting antenna, and to a great extent on 
the nature of the terrain surrounding the transmitter. Similar factors 
affect the energy absorbed from the wave at the receiving station. 

For a station located on or near the ground no energy will be radiated 
below the horizontal, and but little until the angle \f/ of departure of 
the waves above the horizontal becomes appreciable. The minimum 
value of \f/ at which sufficient energy is radiated (or received) to 
produce a readable signal varies with the terrain, the power of the 
transmitter and the sensitivity of the receiver. Over sea water \f/ can 
be very nearly zero; over land the minimum value may be several 
degrees. A fair average approximation may be that ^ must exceed 
about 3.5°. 

A simple geometrical calculation gives yp in terms of sec <£ for various 
distances, and these may be noted on the transmission curves. The 
point of reflection must then correspond to an angle \p greater than 
the minimum value assumed for the terrain at the receiver and 
transmitter. 

The effect of this limitation is to limit the maximum distance for 
single-reflection transmission to about 1,750 km (over land) for the 
£-layer and about 3,500 to 4,000 Ion (over land) for F 2 -l&yer trans- 
mission. Where the transmission path exceeds these values in length, 
calculations must be made on the basis of multireflection transmission. 

The condition of the part of the ionosphere traversed by the wave 
determines the behavior of the wave. In calculating transmission 
conditions, therefore, the ionosphere around the middle of the path, for 
single-reflection transmission, and around the middle of each jump, for 
multireflection transmission, must be considered. The latter require- 
ment means that usually, for practical calculations, one must consider 
the condition of ionosphere over all the transmission path except a few 
hundred kilometers near the transmitter and receiver. This question 
has been discussed in detail elsewhere. 7 It must also be kept in mind 
that at times the o- and ^-components are reflected from different 
geographical parts of the ionosphere. 

VI. DETERMINATION OF SEC O 

Figure 11 gives an alignment chart for the rapid determination of 
the factor sec </> to be used in calculating the logarithmic transmission 
curves. To use the chart, place a straightedge so that it passes 
through the desired virtual height and the desired distance laid off 
on the distance scale at the lower left-hand edge of the chart (in- 
creasing distances lie to the left). The ordinate of the intersection 
of the straightedge with the vertical line corresponding to the same 

7 T. R. Gilliland, S. S. Kirby, N. Smith, and S. E. Reymer. Characteristics of the ionosphere and their 
applications to radio transmission. J. Research NBS 18, 645 (1937) RP1001; Proc. Inst. Radio Engrs. 25, 
823 (1937). 

The weekly radio broadcasts of the National Bureau of Standards on the ionosphere and radio-transmission 
conditions. Letter Circular LC499. 
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desired distance laid off on the main distance scale (increasing 
distances lie to the right) gives the value of sec </> . The 
relation of the point of intersection to the curved dashed lines of 
equal \p gives the value of the angle of departure of the waves from 
the horizontal. A point of intersection falling above the ^=0° line 
indicates an impossible case, where the ray would have to depart 
at an angle below the horizon. 

For example, a distance of 2,400 km and a virtual height of 300 
km corresponds to a sec <j> of 3.07, and an angle of departure of 8.2°. 

VII. TRANSMISSION FACTORS 

When average transmission conditions over a period of time or 
when a variety of transmission paths are to be considered, or when 
an estimate of the maximum usable frequencies is to be made without 
a precise knowledge of the ionosphere over the transmission path, it 
is convenient to have available a means by which the maximum 
usable frequencies may be quickly estimated from an approximate 
value of the vertical-incidence critical frequency. 

The National Bureau of Standards is now beginning a compilation 
of factors by which the critical frequency for the o-component, 
measured at vertical incidence, may be multiplied in order to obtain 
the maximum usable frequencies. These factors are based on average 
observations over a period of time, and may be applied either to 
average critical frequencies to give average transmission conditions, 
or to a given observation of a critical frequency to obtain approximate 
transmission conditions at a given time. 

VIII. CONCLUSIONS 

The type of transmission curves described above have been in 
use at the National Bureau of Standards for the past 2 years in 
studying the correlation of high-frequency radio-transmission con- 
ditions with regular ionosphere observation^. The results of con- 
tinuous measurements of high-frequency broadcasting stations and 
observations on other high-frequency signals, as well as the results 
of some specific experiments have been compared with vertical- 
incidence data. Practically all the available data agree with what 
would be expected on the basis of the theory outlined above, and the 
exceptions may in most cases be accounted for. 

Washington, November 22, 1937. 



